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Abstract 

Nontrivial qq sea effects have their origin in the low-Q 2 dynamics of strong 
QCD. We present here a quark model calculation of the contribution of ss 
pairs arising from a complete set of OZI-allowed strong Y*K* hadronic loops 
to the net spin of the proton, to its charge radius, and to its magnetic moment. 
The calculation is performed in an "unquenched quark model" which has been 
shown to preserve the spectroscopic successes of the naive quark model and 
to respect the OZI rule. We speculate that an extension of the calculation to 
the nonstrange sea will show that most of the "missing spin" of the proton is 
in orbital angular momenta. 
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I. INTRODUCTION 



While providing a good description of low-energy strong interaction phenomena, the 
constituent quark model appears to be inconsistent with many fundamental characteristics 
of QCD. Foremost among these inconsistencies is a "degree of freedom problem": the quark 
model declares that the low energy spectrum of QCD is built from the degrees of freedom of 
spin- 1 fermions confined to a qq or qqq system. Thus, for mesons the quark model predicts 

and we seem to observe a "quarkonium" spectrum. In the baryons it predicts 

and we seem to observe the spectrum of two relative coordinates and three spin-| 

degrees of freedom. 

These quark model degrees of freedom are to be contrasted with the most naive inter- 
pretation of QCD which would lead us to expect a low energy spectrum exhibiting 36 quark 
and antiquark degrees of freedom (3 flavors x 2 spins x 3 colors for particle and antipar- 
ticle), and 16 gluon degrees of freedom (2 spins x 8 colors). Less naive pictures exist, but 
none evade the first major "degree of freedom problem" that the gluonic degrees of freedom 
appear to be missing from the low energy spectrum. This issue, being one of the most 
critical in nonperturbative QCD, is being addressed by many theoretical and experimental 
programs. 

The second major "degree of freedom problem", and the one which we address here, has 
to do with qq pair creation. A priori, one would expect pair creation to be so probable that 
a valence quark model would fail dramatically. Of course, we know empirically that pair 
creation is suppressed: the observed hadronic spectrum is dominated by narrow resonances, 
while the naive picture would predict resonances with widths V comparable to their masses 
m. 

It is now widely appreciated that the narrow resonance approximation can be rationalized 
in QCD within the 1/N C expansion in the limit N c — > oo, meson widths (for example) are 
proportional to iV" 1 while their masses are independent of N c . The demonstration proceeds 
by showing that hadron two-point functions are dominated by graphs in which the valence 
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quark lines propagate from their point of creation to their point of annihilation without 
additional quark loops. The dominance of such a "quenched approximation" is, however, 
not sufficient to underwrite the valence quark model: in the chiral limit, such Feynman 
graphs in general receive important contributions from not only forward quark propagation, 
but also from "Z-graphs" . (A "Z-graph" is one in which the interactions first produce a pair 
and then annihilate the antiparticle of the produced pair against the original propagating 
particle). Cutting through a large N c two-point function at a fixed time therefore would in 
general reveal not only the valence quarks but also a large qq sea. The large N c expansion 
also leaves unanswered a more quantitative question. While hadronic widths F are normally 
small compared to hadronic masses m, they are typically comparable to the mass spacings 
between states in the hadronic spectrum. It is thus surprising that the spectroscopy of a 
valence quark model can survive "unquenching" . 

There is another puzzle of hadronic dynamics which is reminiscent of this one: the success 
of the OZI rule 0. A generic OZI- violating amplitude Aqzi can also be shown to vanish 
like 1/N C . However, there are several unsatisfactory features of this "solution" to the OZI 
mixing problem []3| . Consider uo-cf) mixing as an example. This mixing receives a contribution 
from the virtual hadronic loop process uo — > KK — > 0, both steps of which are OZI allowed, 
and each of which scales with N c like T 1//2 ~ N~ l l 2 . The large N c result that this OZI 
violating amplitude behaves like N~ l is thus not peculiar to large N c : it just arises from 
"unitarity" in the sense that the real and imaginary parts of a generic hadronic loop diagram 
will have the same dependence on N c . The usual interpretation of the OZI rule in this case 

that "double hairpin graphs" are dramatically suppressed is untenable in the light 

of these OZI-allowed loop diagrams. They expose the deficiency of the large N c argument 
since Aqzi ~ T << m is not a good representation of the OZI rule. (Continuing to use u-<f> 
mixing as an example, we note that m^ — m^ is numerically comparable to a typical hadronic 
width, so the large N c result would predict an u-cj) mixing angle of order unity in contrast 
to the observed pattern of very weak mixing which implies that Aqzi « T << m.) 

In our recent papers @J^] , we have studied the unquenching of the quark model, address- 
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ing in particular the impact of qq pair creation on quark model spectroscopy and on the 
OZI rule. In this paper we extend our previous work to calculate some effects of the strange 
quark content of the proton induced by strong ss pair creation. Since, as will be described 
in the next section, our model preserves the spectroscopic successes of the quark model and 
is consistent with the OZI rule, it provides a legitimate framework for the study of the qq 
sea. We focus here on the ss sea both because it allows us to avoid complexities associated 
with antisymmetrization with respect to the valence quarks in the nucleons, and because it 
has recently received considerable experimental attention. 

Our goals for this calculation, though ambitious, are limited. In particular, we will 
address here the effects on the strange quark helicity As, the strangeness charge radius R 2 S , 
and the strangeness magnetic moment fi s of a complete sum over the OZI-allowed ss loops 
which contribute to two-point functions {i.e., of processes that correspond at the hadronic 
level to p — > Y*K* -4 Y*' K*' — > p, where Y* and K* represent generic S = — 1 baryons 
and S — +1 mesons, respectively, and where indicates the action of the appropriate 
current). In contrast, we are unable to discuss the effects of pure OZI-forbidden processes 
(i.e., ones that do not proceed through strong OZI-allowed meson loops). These include 
processes in which the ss pair is directly created or annihilated in a color singlet state (e.g., 
p — * p4>jpc p4>'j lP ic> P and p p4>jpc — > p or p — > p(j)jpc — > p, where 4>jpc is an ss 
meson with quantum numbers J PC . The latter two of these processes correspond to pure 
OZI-forbidden vector-meson-dominance-type graphs. As was the case in our earlier studies 
of OZI violation ||, all such disconnected "double hairpin" diagrams are outside of the scope 
of our model: we focus on the naively much larger OZI-allowed loop diagrams. We also do 
not discuss here processes in which strange baryon-meson loops are directly created by the 
probing current. While such "contact graphs" would in general exist, we show below that 
none are required to make the contributions of our strong Y*K* loop graphs to As, R 2 S , or 
ji s gauge invariant. 
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II. UNQUENCHING THE QUARK MODEL: BACKGROUND 



The Introduction describes three puzzles associated with the nature and importance of 
qq pairs in low energy hadron structure: 

1) the origin of the apparent valence structure of hadrons (since even as N c — > oo, 
Z-graphs would produce pairs unless the quarks were heavy), 

2) the apparent absence of unitarity corrections to naive quark model spectroscopy, 
despite one's expectation of mass shifts Am ~ T (where T is a typical hadronic width), and 

3) the systematic suppression of OZI-violating amplitudes A Q zi relative to one's expec- 
tation (from unitarity) that Aozi ~ r. 

In this section we describe the solutions we see to these puzzles. The resulting picture forms 
the context of the new work described in this paper. 

A. The Origin of the Valence Approximation 

As already mentioned, a weak form of the valence approximation seems to emerge from 
the large N c limit in the sense that diagrams in which only valence quark lines propagate 
through hadronic two-point functions dominate as N c — > oo. This dominance does not seem 
to correspond to the usual valence approximation since the Z-graph pieces of such diagrams 
will produce a qq sea. 

Consider, however, the Dirac equation for a single light quark interacting with a static 
color source (or a single light quark confined in a bag) . This equation represents the sum of a 
set of Feynman graphs which also include Z-graphs, but the effects of those graphs is captured 
in the lower components of the single particle Dirac spinor. I.e., such Z-graphs correspond to 
relativistic corrections to the quark model. That such corrections are important in the quark 
model has been known for a long time |§. For us the important point is that while they 
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have quantitative effects on quark model predictions (e.g., they are commonly held to be 
responsible for much of the required reduction of the nonrelativistic quark model prediction 
that qa = 5/3 in neutron beta decay), they do not qualitatively change the single-particle 
nature of the spectrum of the quark of our example, nor would they qualitatively change 
the spectrum of qq or qqq systems. Note that this interpretation is consistent with the fact 
that Z-graph-induced qq pairs do not correspond to the usual partonic definition of the qq 
sea since Z-graphs vanish in the infinite momentum frame. Thus the qq sea of the parton 
model is also associated with the qq loops of unquenched QCD. 



B. The Am «T Problem 

Consider two resonances which are separated by a mass gap Sm in the narrow resonance 
approximation. In general we would expect that departures from the narrow resonance 
approximation, which produce resonance widths T, ought also to produce shifts Am of 
order T. Yet even though a typical hadronic mass spectrum is characterized by mass gaps 
5m of order 500 MeV, and typical hadronic widths are of order 250 MeV, this does not seem 
to happen. 

We have proposed a simple resolution of this puzzle [f| . In the flux tube model of Ref . [|7| , 
the quark potential model arises from an adiabatic approximation to the gluonic degrees of 
freedom embodied in the flux tube. For example, the standard heavy QQ quarkonium 
potential Vqq(t) is the ground state energy E (r) of the gluonic degrees of freedom in the 
presence of the QQ sources at separation r. At short distances where perturbation theory 
applies, the effect of Nj types of light qq pairs is (in lowest order) to shift the coefficient of 
the Coulombic potential from af\Q 2 ) = 33fa( 1 g /A . ) to af f) (Q 2 ) = (33 , 2jV/) 1 g Q2/A ^ r The 
net effect of such pairs is thus to produce a new effective short distance QQ potential. 

Similarly, when pairs bubble up in the flux tube (i.e., when the flux tube breaks to 
create a Qq plus qQ system and then "heals" back to QQ), their net effect is to cause a 
shift AEjsfJr) in the ground state gluonic energy which in turn produces a new long-range 
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effective QQ potential ||. 

In Ref. Q we showed that the net long-distance effect of the bubbles is to create a new 
string tension b N (i.e., that the potential remains linear). Since this string tension is to be 
associated with the observed string tension, after renormalization pair creation has no effect 
on the long-distance structure of the quark model in the adiabatic approximation. Thus the 
net effect of mass shifts from pair creation is much smaller than one would naively expect 
from the typical width T: such shifts can only arise from nonadiabatic effects. For heavy 
quarkonium, these shifts can in turn be associated with states which are strongly coupled 
to nearby thresholds. 

We should emphasize that it was necessary to sum over very large towers of Qq plus qQ 
intermediate states to see that the spectrum was only weakly perturbed (after unquenching 
and renormalization). In particular, we found that no simple truncation of the set of meson 
loops can reproduce such results. 

C. The Survival of the OZI Rule 

The Introduction illustrates, via the example of u-<f) mixing through a KK loop, why 
unquenching the quark model endangers the naive quark model's agreement with the OZI 
rule. In Refs. |J we showed how this disaster is naturally averted in the flux tube model 
through a "miraculous" set of cancellations between mesonic loop diagrams consisting of 
apparently unrelated sets of mesons (e.g., the KK, KK* + K*K, and K*K* loops tend to 
strongly cancel against loops containing a K or K* plus one of the four strange mesons of 
the L = 1 meson nonets). 

Of course the "miracle" occurs for a good reason. In the flux tube model, where pair 
creation occurs in the 3 Po state, the overlapping double hairpin graphs which correspond 
to OZI- violating loop diagrams (see Fig. 1), cannot contribute in a closure-plus-spectator 
approximation since the ++ quantum numbers of the produced (or annihilated) pair do not 
match those of the initial and final state for any established nonet. Ref. || demonstrates 
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that this approximation gives zero OZI violation in all but the (still obscure) ++ nonet, 
and shows that corrections to the closure-plus-spectator approximation are small, so that 
the observed hierarchy Aozi << T is reproduced. 

We emphasize once again that such cancellations require the summation of a very large set 
of meson loop diagrams with cancellations between apparently unrelated sets of intermediate 
states. 



d d 



(a) 



(b) 



d d 



u u u u 

FIG. 1. (a) OZI-violation in a meson propagator by "pure annihilation", corresponding to a 
disconnected double-hairpin diagram, (b) A different time ordering of the same Feynman graph 
gives an OZI-violating loop diagram via two OZI-allowed amplitudes. 



D. Some Comments 



We believe the preceding discussion strongly suggests that models which have not ad- 
dressed the effects of unquenching on spectroscopy and the OZI rule should be viewed very 
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skeptically as models of the effects of the qq sea on hadron structure: we have shown that 
large towers of mesonic loops are required to understand how quarkonium spectroscopy and 
the OZI rule survive once strong pair creation is turned on. In particular, while pion and 
kaon loops (which tend to break the closure approximation due to their exceptional masses) 
have a special role to play, they cannot be expected to provide a reliable guide to the physics 
of qq pairs. 



III. STRANGE QUARKS AND THE SPIN CRISIS: SOME HISTORY 

Beginning in 1988 with the European Muon Collaboration (EMC) experiment ||, and 
continuing through a recent series of closely related experiments ||10|| , the helicity structure 
functions of quarks in the proton and neutron have been measured via polarized deep inelas- 
tic scattering. When combined with measurements of axial charges in hyperon beta-decay 
and the assumption of SU(3) symmetry, these experiments indicate a "spin crisis": only 
about a third of the nucleon's helicity resides on its quarks, and about —10 ± 3% of this 
helicity is lost to strange quarks in violation of the Ellis- Jaffe extension \Y1\ of the 
fundamental Bjorken Sum Rule [|T3j]. 

Although generally accepted, there has been some discussion about the reliability of 
these conclusions. While support for them has come from other types of experiments ||14|| , 



they have been criticized from other quarters for depending on an extrapolation of the 
structure functions to small x [15| and on an S'L r (3)-symmetry-based analysis of hyperon 



beta decay [[T(J. At a deeper level, reanalyses of the theoretical connection between spin- 
dependent deep inelastic scattering and the spin structure functions showed that the SU (3) 
singlet structure functions are entangled with the gluon spin structure functions via the 
£7(1) axial anomaly |[7). This observation has led to attempts to avert the "spin crisis" by 
invoking a large gluonic contribution via the anomaly. This possibility should be checked 
by direct measurements on the glue. 

The naive nonrelativistic quark model predicts that 100% of the nucleon's helicity resides 
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on the quarks, but, as already mentioned above, lower components of the quark spinors 
arising from relativistic effects are believed to lower this fraction to about 75% H]. At the 



opposite extreme are naive Skyrmion models [|18|] which predict that the net quark spin of 
the nucleons should be zero (a result which seemed supported by the initial experimental 
results) . 

If there is a large strange quark contribution to the nucleon spin, then one would also 
naturally expect strange contributions to nucleon magnetic and electric form factors. Purely 
electromagnetic scattering can only measure the four linear combinations 

n yp _ 2 r y(u) 1 Md) _ ^ Ms) /.x 
^M,E — '^ jr M,E ~^ M > E M,E \ L ) 

9 11 

( f) 

where G^J E is the magnetic (M) or electric (E) form factor of the quark flavor / in the 
proton. From parity violating scattering on the proton one can measure two more linear 
combinations 

GHe = (\- \sin*6 w )G% + (-i + \smH w )[G^ E + G%] (3) 

and thereby separate out the six elementary form factors G^ E for / = u, a, and s. Experi- 
ments are currently underway and others are planned to measure these form factors. Such 
measurements appear to be the next step in understanding the physics of the spin crisis. 

In the wake of the spin crisis have come a number of attempts to find theoretical descrip- 
tions less extreme than the naive quark and Skyrmion models. In 1989, Jaffe [19| pointed out 



that the pole fit of Hohler et al. to the nucleon's isoscalar electromagnetic form factors f20 



could suggest the presence of significant strange currents in the nucleon. By identifying the 
two lightest fitted poles with the physical oj and cj) mesons, he estimated R 2 S = 0.14±0.07 fm 2 
for the strangeness radius and fi s = —(0.31 ± 0.09) /j,n for the strange magnetic moment of 
the nucleon. (Note that Jaffe uses a sign convention in which the strange quark has positive 



strangeness, opposite to the Particle Data Group's convention | 21| . As a result, both R 2 S 
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and fi s must be multiplied by — | to obtain the contributions of the strange quarks to the 
electric and magnetic form factors of the proton. Thus, for example, a positive value for 
R 2 S indicates that s quarks are farther on average from the proton's center than s quarks. 
Jaffe's convention appears to have been adopted by subsequent authors, and we too shall 
adhere to it.) 

Jaffe and Lipkin [p2], building on earlier work by Lipkin p3fl , constructed an extended 



quark model in which the valence qqq component of octet baryons was supplemented with a 
"sea" consisting of a single qq state which was allowed to have either ++ or 1 ++ quantum 
numbers. Their model was not predictive; it was intended only as an example of a simple 
extension to the quark model which could accommodate the EMC results as well as baryon 
magnetic moments and hyperon /3-decay. They found that the data could be fit with either 
a {uu + dd + ss) or {uu + dd) flavour structure to the sea, though in both cases a large 
suppression of the purely valence component of the baryon wavefunctions was required. For 
other early analyses along these lines, see Refs. |24| . 

The renormalization of axial couplings qa (and therefore of the fraction of the proton 
spin Aq carried by the quarks of flavor q) by qq pairs in the form of meson loops is a subject 
with a history dating back to the birth of meson exchange theories of the strong interaction. 
For some modern studies in the context of chiral perturbation theory, see Ref. [pql . Many 
recent studies, including ours, are extensions of this classic meson loop approach ||26|| . 

A model-dependent study of the ss sea based on hyperon-kaon loop diagrams was made 
by Koepf et al. in Ref. ||27|| . These authors used both a non-relativistic quark model and 
the cloudy bag model to calculate the strangeness content of the proton arising from AK, 
UK, and Y?K loops. After tuning the baryon-baryon-meson form factors to reproduce the 
nonstrange nucleon moments, they found that both models predict rather small strange 
moments: As ~ —0.003, R 2 S ~ — .Olfm 2 , and fi s ~ — .03fi N . 

Subsequently, Musolf and Burkhardt EH] examined the AK loop graph in a calculation 



which took its vertex form factors from the Bonn potential for baryon-baryon scattering [R0 



and which included seagull graphs in order to satisfy the vector current Ward-Takahashi 
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identity. The results obtained by these authors, As ~ —0.044, R 2 ~ — .03fm 2 , and /x 



— 0.35/itv, are significantly larger than those found in p7| . The discrepancy is due, at least 
in part, to the non-gauge- invariance of the earlier calculation. 

More extensive, but from our perspective still incomplete, meson loop calculations have 
been presented in Refs. [[H]]. These authors extend loop calculations to include the entire 
set of ground state octet pseudoscalar and vector mesons plus the ground state octet and 



decuplet baryons. They obtain As ~ +0.02. Another tack has been taken by Ito ||29|| , who 
calculates the effects of kaon loops on a constituent quark, obtaining R 2 w — .02fm 2 and 
fi 8 « -.12 n N . 

The loop and pole pictures were combined in the model of Cohen et al. [32| , who obtained 
R 2 ss -.042 fm 2 and Us ~ —.28 hn- These authors also calculated the strange form factors 
at nonzero momentum transfer. 

For some reviews and for alternative models and points of view, in particular Skyrme- 



based calculations, see Refs. 33 



IV. A PAIR CREATION MODEL FOR THE STRANGENESS CONTENT OF 

THE PROTON 

Our discussion of the strangeness content of the proton will be based on the quark-level 
process shown in Fig. ||](b). The main new feature of our calculation is that we shall sum 
over a complete set of strange intermediate states, rather than just a few low- lying states. 
Not only does this have a significant impact on the numerical results for As, R 2 , and /i s , 
but, as explained above, it is necessary for consistency with the OZI rule and the success of 
quark model spectroscopy. 
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FIG. 2. A meson loop correction to a baryon propagator, drawn at (a) the hadronic level, and 
(b) the quark level. 



The lower vertex in Fig. 0(b) arises when qq pair creation perturbs the initial nucleon 
state vector so that, to leading order in pair creation, 

|P>-|P>+ E h 2 dd \Y*K*qlS) { l* K J Slh ^ P) , (4) 

Y*K*IS LVLp — — &K* 

where h q q is a quark pair creation operator, Y* (K*) is the intermediate baryon (meson), q 
and I are the relative radial momentum and orbital angular momentum of Y* and K*, and 
S is the sum of their spins. Of particular interest is ss pair creation by the pair creation 
operator h S g, which will generate non-zero expectation values for strangeness observables: 



Y*'K*'e'S' 

(5) 

The derivation of this simple equation, including the demonstration that it is gauge invariant, 
is given in the Appendix. We will be considering the cases O s = As, R 2 S) and /x s . The value 
of As can be associated (via small scale-dependent QCD radiative corrections) with the 
contribution of strange quarks to the deep inelastic spin-dependent structure functions and 
to the strange quark axial current matrix elements in the proton. 

To calculate the p — > Y*K* vertices in Eq. (Q), we employ the same flux-tube-breaking 
model used in our earlier work. This model, which reduces to the well-known 3 P decay 
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model in a well-defined limit, had its origin in applications to decays of mesons [0,0 and 
baryons ||36|| . The model assumes that a meson or baryon decays when a chromoelectric flux 
tube breaks, creating a constituent quark and antiquark on the newly exposed flux tube 
ends. The pair creation operator is taken to have 3 Pq quantum numbers: 



3/2 



hqq{t, X) 



7o 



J d 3 z exp 



3z 2 



g t (t,x + -)a- Vg(t,x--) . 



(6) 



The dimensionless constant 70 is the intrinsic pair creation strength, a parameter which 
must be fit to decay data. In our previous studies of mesons, we fit 70 to the p — > tttt decay 
width; here we find it more appropriate to fit to the A — > Ntt width. It turns out that 
the two values agree to within 20%, which is a reassuring consistency check. The operator 
(§) creates constituent quarks, hence the pair creation point is smeared out by a gaussian 
factor whose width, r q , is another parameter of the model. (In addition to being physically 
motivated, this smearing factor is required to render the sum in Eq. (|5]) finite.) As discussed 
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r„ is constrained by meson decay data to be approximately 0.25 fm. 



Once an ss pair is created, the decay proceeds by quark rearrangement, as shown in 
Fig. [| The p — > Y*K* decay amplitude of the first of Figs. ^ may be written as 



(Y*K*\h s - s \p)= l0 X-I, 



(7) 



where £ is a spin overlap which can be expressed in terms of the baryon and meson spin 
wavef unctions as 



(8) 



with 



XS4S5 



-25 Sil 5 S5l j 



(9) 



and / a spatial overlap: 
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p 



m,. 



-q 



pexp ( ~-r 2 q p 2 ) $ p 



2V P+ 6" q 



(10) 



Here the $'s are momentum space wavef unctions, q is the momentum of Y*, and the m/s 
are quark masses {m uus is short for 2m u + m s , etc.). The factor exp(— s 2 /2b) models the 
overlap of the initial and final-state flux tube wavef unctions; its size is controlled by the 
physical string tension b, though our results do not depend strongly on the numerical value 
of b. 

For the remaining quark line diagrams in Fig. |3], the decay amplitude still has the form 
(0), but the spin indices in Eq. (D) become permuted. (The spatial overlap in (|T0|) remains 
the same thanks to the assumed symmetry of the proton's spatial wavef unction.) 



u u 



u u 



u d s 



1 1 



u d s 



1 1 







FIG. 3. Quark line diagrams for p — ► T,*K* and p — > A*K* 



Faced with the large number of states that contribute to the sum in Eq. ([5]), we have 
found it necessary to use simple harmonic oscillator (SHO) wavefunctions for the baryons 
and mesons in (P]). The oscillator parameters (3 (defined by $(k) ~ e~ fc2 / 2/32 ), were taken to 
be Pmeson = 0.4 GeV for mesons (as in Ref. H35[] ) and Pbaryon = 0.32 GeV for baryons (as in 
Ref. f37|). As discussed below, our results are quite insensitive to changes in the /3's (mainly 

because Eq. ([5D is independent of the choice of wavefunctions in the closure limit any 

complete set gives the same result and the full calculation with energy denominators 
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does not deviate much from this limit.) 

Even with SHO wavef unctions, the sum over intermediate states would be very difficult 
were it not for an important selection rule: inspection of the quark line diagrams in Fig. |3| 
shows that the relative coordinate of the non-strange quarks in baryon Y*K* is always in 
its ground state. Only the relative coordinate between the strange and non-strange quarks 
{i.e., the Ay. -oscillator) can become excited. This drastically reduces the number of states 
that must be summed over. Unfortunately, this simplification does not apply for uu or dd 
pair creation; we therefore postpone to a later paper the computation of their contributions 
to the nucleons' spin, charge radii, and magnetic moments. 

We will find it useful at times to refer to the closure-spectator limit of Eq. (||). This is 
the limit in which the energy denominators do not depend strongly on the quantum numbers 
of Y* and K*, so that the sums over intermediate states collapse to 1, giving 

(O s ) oc (p\h sS O s h sS \p) oc (0 \h S sO s h sS \ 0) , (11) 

where the second step follows since h S s does not couple to the motion of the valence spectator 
quarks. We see that the expectation value of O s is taken between the 3 Pq states created by 
h sS . From the properties of the 3 Pq wavefunction it then follows that As = R 2 = fi s = 
in the closure-spectator limit (a result which would not be seen if only the lowest term, or 
lowest few terms, were included in the closure sum). 

In the next Section we will present our results for the expectation values defined by 
Eq. (|5|) for the quantities As, R 2 , and We will see that delicate cancellations lead to 
small values for these observables even though the probability of ss pairs in the proton is of 
order unity! 

V. RESULTS 
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A. Strange spin content 



As, the fraction of the proton's spin carried by strange quarks, is given by twice the 
expectation value of the s and s spins : 

As = 2 (Sl s) + S { z s) ) . (12) 

Let us first examine the contribution to As from just the lowest-lying intermediate state, 
AK. The P-wave AK state with J = J z 



(AK), 



m=l 



(A T JO 



m=0 



(13) 



The s quark in the kaon is unpolarized, while the s quark in the A carries all of the A's spin; 
because of the larger coefficient multiplying the first term in (|T3"D, the AK intermediate state 
alone gives a negative contribution to As. 

When we add in the (AK*) P i and (AK*) P 3 states (note that the subscripts denote the 
quantities £S defined previously), we have 



As oc (1 -J\ 



1 



-3 

-1 



10 



/ 1 \ 



V 



(14) 



in the closure limit. Here the matrix is just 2(S^ + S^') (which is of course symmetric 
though we only show its upper triangle), and the vectors give the relative coupling strengths 
of the proton to [ (AK) P i, (AK*) P i, (AK*) P z }. There are a couple of things to note here: 
(1) The matrix multiplication in (|l"4D evaluates to zero; there is no net contribution to 
As from the AK and AK* states in the closure limit. There are in fact many such "sub- 
cancellations" in the closure sum for As: for each fixed set of spatial quantum numbers in 
the intermediate state, the sum over quark spins alone gives zero (because (S^) = (S^) = 
in the 3 Po state). That is, each SU(6) multiplet inserted into Eq. (|5|) separately sums to 
zero. Moreover, the As operator does not cause transitions between 1 = and 1 = 1 strange 
baryons so that the A and E sectors are decoupled, hence they individually sum to zero. 
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(2) Only the diagonal term in Eq. (|TJJ) corresponding to p — > (AK*) P s ^(AK*) P 3 — > p 
gives a positive contribution to As. (We use — * here to denote the action of the As opera- 
tor.) All of the other terms give negative contributions. In the full calculation with energy 
denominators, the negative terms are enhanced because they contain kaon (rather than K*) 
masses. The full calculation gives As = —0.065 from AK and AK* states. The largest 
individual contribution is —0.086, from the off-diagonal term p — > [AK) P i— >(AK*) P 3 — > p. 
Proceeding to intermediate states containing E and E* baryons, we calculate 

^3 -12^2 3^ -6^6 
6^ 

-7 -10V2 
10 



2(S« + S<?>) 



1 

54 



15 



(15) 
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in the basis [ (EK) p i, (E*fC) P 3, (Eif*) p i, (Eif*) p |, (E*if*) p i, (E*i^*) p 2 ]. The corre- 
sponding relative couplings to the proton are 

Again, the net As from these states is zero in the closure limit, but this time the insertion 
of energy denominators does not spoil the cancellation very much: the full calculation gives 
As = —0.003 in this sector. 

P-wave hyperons and kaons contribute another —0.04 to As, and the net contribution 
from all higher states is —0.025. Thus, the result of our calculation is As = —0.13, in quite 
good agreement with the most recent extractions from experiment, As (expt) = — 0.10±0.03 
(see, e.g., Ref. fill). We emphasize that our parameters were fixed by spectra and decay 
data. Moreover, the result is quite stable to parameter changes, varying by at most ±0.025 
when r q , b, /3b ary0 n and Aneson are individually varied by 30%. 

For comparison with other calculations, we note that in our model the AK intermediate 
state alone contributes —0.030 to As, and the contribution from the AK, UK, and H*K 
states together is (coincidentally) also —0.030. 

It is interesting to note that As is driven mainly by meson, rather than baryon mass 



splittings: if we set m\ = = m^., we find that As decreases by only about 30%, whereas 
it drops by about 80% if we set = m^*. Finally, we have also computed the charm-quark 
contribution to the proton spin, finding Ac ~ —0.01. 



B. Strangeness radius 

Figure § defines our variables for the quarks in an intermediate state. The (squared) 
distances of the s and s quarks from the baryon-meson center of mass are 



(r 4 - R c 



-VQ ( ) Ay. + e K *r 

\m„, 



(16) 



(r 5 — R cm ) 2 



m,. 



v K * - e Y *r 



(17) 



where €k* = Mk* / 'My*k* and ey* = My* / 'My*k*, while by definition 



Rt = rt- r 



(18) 



is the strangeness radius. 




FIG. 4. Quark coordinates in an intermediate state. 
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The calculation of R 2 is more difficult than the calculation of As, for several reasons. 
First, the operators appearing in R 2 cause orbital and radial transitions among the interme- 
diate states. Thus SHO transitions satisfying An = 0, ±1 and/or A£ = 0, ±1 are allowed, 
so there are many more terms to calculate (n and £ are orbital and radial SHO quantum 
numbers). Moreover, the sub-cancellations discussed above no longer occur, so that R 2 S con- 
verges more slowly than As: we must include more states in Eq.(f|) to obtain good accuracy. 
In addition, the basic matrix elements are more complicated: in a basis of states with good 
magnetic quantum numbers (m, m') we have, for example, 



(n'£'m'\rK*z\n£m) 



(£ + m)(£-rn) 



\ (2£+l)(2£ 



(jn + £ + l/28 n , n - Vn + 15 n , n+1 ) 



(£ + m + l)(e-m+l) ( I ; ~ - . ^ r \ 

\ (2<+l)(2< + 3) (V» + < + 3/^ " " ^ »-.) 



(19) 



for matrix elements of the meson internal coordinate and 



(q'£'m'\r z \q£m) = i8 m > m \ 8 e e . 



(£ + m)(£ - m) 
\ {2£+ l){2£- 1) [ dq 



d £-1 
— + 



Q 



8(i 



l' £+1' 



(£ + m + l)(£-m+ 1) 
\ (2£+l)(2£ + 3) 



d_ £ + 2 
dq q 



8{q - q') 



(20) 



for matrix elements of the Y* — K* relative coordinate. These matrix elements must be 
coupled together to give {R 2 S ) between states of definite £ and S with total angular momentum 
|, leading to formulas which become quite lengthy, especially for excited intermediate states. 
Thus we were happy to have a stringent check of our results: when we equate all of the energy 
denominators in Eq. @, we must obtain the closure-spectator result, R 2 = 0. 

Our results for R 2 are shown in Table |I[ With our standard parameter set, we obtain R 2 = 
— 0.04fm 2 . For reasonable parameter variations, R 2 ranges between —0.02 and — 0.06fm 2 . 
Table | shows that the lowest lying SU(Q) multiplets of intermediate states (i.e., the S"-wave 
hyperons and kaons) account for about half of r 2 and r 2 . Most of the remaining contributions 
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come from P-wave hyperons and kaons. However, R 2 involves a large cancellation between r 2 
and r 2 , and its value doesn't settle down until we add in quite highly excited intermediate 
states. For this reason, the precise numerical value (and perhaps even the sign) of R 2 
cannot be considered definitive: our conclusion is rather that R 2 is small, about an order of 
magnitude smaller than r 2 and r 2 . This result is not too surprising: R 2 is exactly zero in 
the closure limit, and our previous hadronic loop studies [f||| led us to expect that the full 
calculation with energy denominators would preserve the qualitative features of this limit. 

Note that the AK intermediate state alone gives R 2 S « -O.Olfm 2 (the sign is as expected 
from the usual folklore) while the AK, Y,K, and ~E*K states together give -0.017fm 2 . Nev- 
ertheless, although our sum over all states gives the same sign and order of magnitude as 
these truncations, Table p] shows that this is just a coincidence. 



TABLE I. Proton strangeness radius from hadronic loops (in fm 2 ). The rows give the running 
totals as progressively more excited intermediate states are added into the calculation. The final 
column thus shows the total from all intermediate states. 
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C. Strange magnetic moment 



The strange and antistrange quarks carry magnetic moments — tt/i^' 5 ) where 



2m 



// 



(s) 



(21) 



(22) 



and we denote the net strange magnetic moment by /x s : 



(23) 



The spin expectation values are already in hand from our As calculation. Referring again 
to Fig. £|, we see that the s and s orbital angular momenta are given by 



L s = (r 4 - R cm ) x p 4 
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(25) 



Computing the expectation values of these operators presents no new difficulties beyond 
those encountered in the R 2 S calculation. In fact, there are no radial transitions in this case, 
so there are fewer states to sum over and the sum converges more quickly. 
The results obtained with our standard parameter set are 

(2S^) = -0.058 (2S&) = -0.074 

(LW) = 0.043 (4 s ")) = 0.038 

fi^ = -0.025/i7v = 0.060/i W 

fj, s — 0.035/iAr 

We predict a positive (albeit small) value for /j, s , in disagreement with the other models dis- 
cussed at the beginning of this section. Where does the positive sign originate? First note 
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that the signs of (Sjf'), (L^), and (L^) are correctly given by just the lowest lying interme- 
diate state, AK of Eq. (|H|). (We also note in passing that the L z 's have similar magnitudes so 
that orbital angular momentum contributes very little to fi s in any case.) On the other hand, 
the AK state has (S^) = 0, whereas we find (S^) to be quite large and negative. (The main 
contribution comes from the off-diagonal process p — > (AK) P i -^(AK*) P 3 -> p, although 
there is also a significant contribution from p — > (A(1405)if) 5 i (A (1405) If*) s i — > p.) 
These important terms, which drive fi s positive, are omitted in calculations which include 
only kaon loops. (We find that the AK intermediate state alone contributes — 0.080/iat to 
fi s , and the contribution from AK, T,K, and E*if together is — 0.074//jv) 

VI. CONCLUSIONS 

We have presented here parameter-free calculations of the effects of the ss sea generated 
by strong Y*K* loops on the low energy, nonperturbative structure of the nucleons. These 
calculations represent what is to our knowledge the first such results within a framework 
which has been demonstrated to be consistent with the many empirical constraints which 
should be applied to such calculations, namely consistency with the success of the quark 
potential model and especially with the validity of the OZI rule. 

Our results predict that observable effects from the strange sea generated by such loops 
arise from delicate cancellations between large contributions involving a suprisingly massive 
tower of virtual meson-baryon intermediate states. If correct, our conclusions rule out the 
utility of a search for a simple but predictive low energy hadronic truncation of QCD. While 
complete (in the sense of summing over all OZI-allowed Y*K* loops) and gauge invariant, 
we recall that our calculation has ignored pure OZI-forbidden effects as well as those loop 
diagrams directly generated by the probing current (contact terms). As a consequence, our 
results cannot strictly speaking be taken as predictions for As, R 2 s1 or fi s . Rather, this 
calculation shows that a complete set of strong Y*K* loops, computed in a model consistent 
with the OZI rule, gives very small observable ss effects. While such OZI-allowed processes 
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might dominate, we cannot rule out the possibility (as was also the case with u — <fi and 
other meson mixing ||) that direct OZI violation (and in this case contact graphs as well) 
could make additional contributions of a comparable magnitude. 

The small residual effect we calculate for the loop contributions to As seems consistent 
with the most recent analyses of polarized deep inelastic scattering data. Our calculations 
also give small residual strange quark contributions to the charge and magnetization distri- 
butions inside the nucleons. If these contributions are dominant, it will be a challenge to 
devise experiments that are capable of seeing them. Indeed, they are sufficiently small that 
we would expect that their observation will require the development of special apparatus 
dedicated to this task. Given the fundamental nature of the puzzling absence of other signals 
for the strong qq sea in low energy phenomena, this effort seems very worthwhile. 

It would be desirable to devise tests of the mechanisms underlying the delicate cancel- 
lations which conspire to hide the effects of the sea in the picture presented here. It also 
seems very worthwhile to extend this calculation to uu and dd loops. Such an extension 
could reveal the origin of the observed violations [38] of the Gottfried Sum Rule [39| and 
also complete our understanding of the origin of the spin crisis. From our previous calcu- 
lations [0, the effects of "unquenching" strange quarks are a good guide to the effects to 
be expected from up and down quarks in the absence of Pauli blocking. Since most of the 
created pairs are in highly excited states, Pauli blocking should be of minor importance, 
and so one would guess that each of up, down, and strange will produce a contribution to 
Aq of about —0.13. When combined with the relativistic quenching of Aq va i ence 0, this 
makes it plausible to us that most of the "missing spin" of the proton is in orbital angular 
momentum. 
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APPENDIX: DERIVATION OF (O s ) 

While Eq. (5) has a simple nonrelativistic interpretation as an expectation value of O s in 
the dressed state of Eq. (4), is it not obvious that it actually computes As, R 2 , and ji s . In 
particular, R 2 and /i s are defined in terms of charge and magnetic form factors, respectively. 
Here we show that the contributions to all three quantities from ss pairs arising from strong 
Y*K* loops may indeed be calculated in a gauge invariant fashion in our model via this 
simple formula, with each given by the usual nonrelativistic operators of Eqs. (12), (18), 
and (23). 

We first discuss the relatively simple case of As, where gauge invariance plays no role. To 
second order in h sS , the graphs contributing to the nucleon matrix elements of = S7^7 5 s 
are of two types: 1) OZI-allowed graphs where a strong ss loop has created a Y*K* loop 
and in which scatters the s or s quark, and 2) pure OZI-forbidden graphs in which 
A^ acts on a color singlet ss state which is created and/or destroyed on the nucleon. The 
pure OZI-forbidden processes fall outside the scope of our model. On the other hand, OZI- 
allowed contributions may be calculated (with all the usual caveats and approximations of 
nonrelativistic quark model calculations of axial matrix elements) according to 

g ^ = ( p (0,+)\7sYl 5 s\p(0,+)) (27) 
^2(p(0,+)|5<') + 5W|p(0,+)> (28) 
= As (29) 

as used in Section V.A. 

The discussion of R 2 S and /i s is considerably more involved. Our first step is to set up 
the Breit-frame formalism for calculating R 2 S and /i s . It is easily shown that the electric and 
magnetic form factors of the proton may be calculated via the formulas 



= -Q 2 ) = (») 
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Q 2 ) 



e<(Q) 



(31) 



where T|" 



(Q) are the T-matrix elements 



n xt (Q)^{p(+^+)\T E \ P ( 




(32) 



T e M xt (Q) = ( P (+^,-)\TM\p( 




(33) 



for scattering in the external potentials 



/ 



rf 3 g e tQ - x 0^(Q) 



(34) 



and 




if) J d 3 Qe^ x <(Q) 



(35) 



respectively. (Note that all normalizations are fixed by the condition that G P E (0) = 1.) 

Now consider a proton scattering off one of these external potentials subject to the 
additional gauge invariant perturbation created by ss pair creation. The resulting changes 
AG P E and AG P M are and G$, which may therefore be associated with ATf^ and AT^f*. 
Thus, generically, the AG's may be associated with four processes. Associated with the 
"naive" formula (5) are the two processes in which: 

la) the external field probes the K* of the Y*K* loop, and 
lb) the external field probes the Y* of the Y*K* loop. 

We will recall at the end of this Appendix how these processes, which are of second order in 
the strong ss pair creation Hamiltonian density h sS and third order in perturbation theory, 
lead to Eq. (5). 

Associated with h sS , which has matrix elements proportional to (— l) 2s xtvX-s- (p s ~ Ps)> 
is a pair creation contact interaction density h c s ° s ntact with matrix elements proportional to 
2e s (— l) 2s xJo"X-s • A = j c s ° ntact ■ A, where e s = — |e and A is the vector potential. Through 
this contact interaction two more processes can occur: 
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2a) the external field directly creates the ss pair through h c ™ tact , converting the proton 
at momentum — ^p into a Y*K* system at momentum +^p, and then /i s5 annihilates the ss 
pair, returning the system to a proton at the same momentum +^p, and 

2b) converts the proton at momentum — ^p into a system at the same momen- 

tum, and then h c s ° s ntact annihilates the ss pair, turning the Y*K* system at momentum — ^p 
to a proton at momentum +^p. 

We begin by showing that the contact terms 2a) and 2b) are not required to make the loop 
contributions la) and lb) to B? s or /i s gauge invariant. 

The case of R 2 S is trivial: h c s ° s ntact is independent of 0f?*(t,x) in our model, so contact 
processes do not contribute to it. 

In contrast, ji s is determined by scattering in A^'(t, x), which does couple to h c jp tact . 
From ordinary second-order time-ordered perturbation theory 



2Tjg*(Q) 
ea^(Q) 




where A(Q 2 ) is an energy denominator which can easily be shown to be even in Q, and 
where 

tff = ^(o,o)|[y*K*]^(+^,-))([y*K*]^(+^,-)|j— c *(o,o)_ (37) 

and 

tf = jr rt (o 1 o)-|[rin^(-f ) +))([rrw-| ) +)|Mo ) o) , (38) 

where \\Y*K*] q es(±Qf,TnY*K*)) denotes a state with internal radial momentum q, internal 
orbital angular momentum £ and total internal spin S coupled to total angular momentum 
and parity J p = | + (so that h sS can connect it to the proton), with momentum ±^p, and 
with ^-component of spin my*K*- Now consider the general matrix element 

([Y*K% es (+^,m Y , K ,)\h s ,(0,0)\p(-^,m p )) = g q Y i S K * p (Q 2 )x Y *K*X P • (39) 

Since h sS is a scalar operator, the form of the right-hand side of this equation is uniquely 
determined. In particular, it follows from parity conservation that this matrix element has no 
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spin-flip component and therefore that when made gauge invariant via contact interactions, 
it will not generate one. We conclude that also no contributions to /i s from contact terms 
are required to make the contributions of the strong Y*K* loops gauge invariant. 

There is, of course, no reason why the underlying ss pair creation dynamics might not 
generate independent gauge invariant contact currents of the form XY*K* e ^ k(7 ^Q k Xp which 
could contribute to /i s . Although they lie outside our goal of providing a complete calculation 
of the effects of ss pairs generated by strong Y*K* loops, it would be interesting to know the 
size of such effects. (We note in passing that, as a class of diagrams, they will also be subject 
to strong cancellations since, for example, 3 Si pair creation followed by 3 Po annihilation will 
vanish in the closure-spectator approximation). 

To complete our demonstration that Eq. (5) is the correct gauge invariant formula for 
the contributions of ss pairs from strong Y*K* loops to As, R 2 S , and /i s , we briefly recall 
how the third-order processes la) and lb) lead to it. The analysis is very straightforward, 
giving 

T e M xt (Q) la)+lb) = T e M xt (Q) Y * + T™\Q) K * (40) 

corresponding to scattering from intermediate particles Y* and K*, respectively, with, for 
example, 

T e M xt (Q) K * = ej£*(Q)<(Q) £ / d\¥£ Y ; mK ^)^ mK ^ + Cy .Q) (41) 

m Y *m K * 

= e^(Q)<(Q) £ /rf'^.-W^C.+W ( 42 ) 

m Y * m K * 

where r^* = ey*r and where 

= ([Y\e Y ,P + n,m Y ,)K%e K ,P -n,m K ,)] es \h s ,(0,0)\p(F,s)) (43) 

is the internal relative momentum wavefunction of Y* and K* "inside" the proton (which 
is independent of the (small) total momentum P), and ip is its Fourier transform. Clearly 
an exactly analogous formula appears for T|f*(Q)y* and the two, having as they do simple 
constituent interpretations, lead to the "naive" Eq. (5). 
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